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ABSTRACT 


This  report  represents  an  extension  of  the  work  described  in  an 
earlier  report  [Electromagnetic  Fields  in  the  Ocean  Near  a  Shoreline, 

J.E.  Spence,  E*  Sullivan  and  J.  Beville  -  Contract  Nonr  396(10)- 
Report  396(lO)/2  -  15  November  1965].  It  is  an  analytical  study  of 
the  electromagnetic  field  in  the  ocean  resulting  from  excitation  over 
land  by  natural  electromagnetic  noise  such  as  tnicropulsations  or  ELF 
atmospherics.  The  main  purpose  of  the  study  is  to  determine  the  rela¬ 
tive  importance  of  electromagnetic  energy  entering  the  ocean  via  the 
air-water  interface  versus  energy  entering  via  the  soil  and  ocean  bottom. 
Emphasis  in  this  report  is  placed  on  a  comparison  between  the  vertical 
and  horizontal  components  of  the  electric  field.  Vertical  profiles  of 
the  horizontal  and  vertical  electric  field  vectors  are  plotted  for 
several  values  of  frequency,  ocean  depth,  and  the  distance  from  the  shore¬ 
line.  Finally,  a  comparison  is  made  with  some  experimentally  determined 
values  of  the  horizontal  electric  field. 
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I.  INTRODUCTION 


Recent  experitaortal  measure-.ner.ts  of  electromagnetic  noise  in  the  ocean 

have  sometimes  produced  data  which  show  an  unexpected  intensity  distribution 

[123 

with  depth  near  a  coastline  .  In  this  thesis  the  electromagnetic  field 
which  may  be  produced  in  the  ocean  by  natural  sources  over  land  is  investi- 

t 

gated  analytically. 

A  simplified,  two-dimensional  model  is  used  which  assumes  a  constant 
depth  ocean  with  an  abrupt  shoreline.  The  field  in  the  ocean  is  assumed 
to  be  excited  by  electromagnetic  energy  passing  across  the  coastline,  into 
the  s ea.  The  magnetic  vector  is  taken  to  be  parallel  to  the  shoreline. 

A  superposition  method  is  used  which  assumes  the  field  in  the  sea 
to  be  composed  of  two  fields  -  one  produced  by  energy  passing  through  the 
air-ocean  interface  and  the  other  produced  by  energy  passing  through  the 
•oil  and  ocean  bottom.  The  energy  passing  through  the  vertical  land-ocean 
interface  ia  neglected  in  thla  treatment. 

The  waves  passing  through  the  air- oca an  interface  art  assumed  to  be 
produced  by  energy  which  passes  across  the  coastline  on  the  surfaces  of  the 
earth  and  the  ocean.  The  ocean  ia  first  assumed  to  have  infinite  conductivity 

which  allow?  the  waves  to  enter  the  ocesn  with  e  direction  of  propagation 

»  • 

normal  to  the  sea  surface.  The  assumption  of  infinite  conductivity  la 
then  dropped.  The  energy  passing  through  the  earth-ocean  interface  Is  treated 
in  a  somewhat  similar  manner,  the  difference  being  that  the  fiald  Intensity 
at  the  ocean  bottom  ta  determined  by  an  application  of  Creen's  integral  theorem, 

! 


I 


The  use  of  Green's  theorem  results,  as  is  usually  the  case,  in  an  in¬ 
tegral  equation  which  must  be  solved  by  an  approximate  method.  In  the  pre¬ 
sent  case  the  solution  leads  to  an  integral  which  is  solved  by  a  far-fleld 
approximation.  The  resulting  expressions  for  the  electric  field  components 
are  examined  through  the  use  of  inequalities  and  the  vertical  component 
of  the  electric  field  vector  is  found  to  bo  negligible  in  comparison  to 
the  horizontal  component  under  the  restrictions  of  this  model. 

The  horizontal  component  of  the  electric  field  is  plotted  for  several 
values  of  frequency,  depth  of  ocean,  and  distance  from  the  shoreline.  A 
comparison  is  then  made  with  some  experimentally  determined  values  of  the 
horizontal  electric  field. 


II.  FORMULATION  OF  THE  PROBLEM 


«**• 

The  first  pert  of  this  chapter  is  concerned  with  a  description  of  the 
ocean*- shoreline  model.  The  forr.;  of  the  electromagnetic  field  incident  on 
the  shoreline  is  then  presented.  Finally,  a  mathematical  expression  is 
obtained  which  describes  the  magnetic  field  in  the  ocean  due  to  the  inci¬ 
dent  field.  i 


2.1  Ocean-Shoreline  Model 


The  geometry  of  the  simplified  two-dimensional  ocean- shoreline  model 
is  shown  in  Figure  1.  In  this  model  the  incident  field  is  assumed  to  have 
a  direction  of  propagation  in  the  ♦  x  direction  and  all  fields  are  con¬ 
sidered  to  be  independent  of  a,  i.et,d/dz  •  0. 

2.2  Incident  Field 

The  magnetic  field  intensity  H  will  be  considered  to  have  a  z  com- 
ponent  only,  Vith  this  restriction,  and  the  assumption  that  d/i s  ■  0, 
Maxwell's  curl  H  equation  yieldsl 


B  • 


ah 


x  cr-iotf  Ey 


.  -i  ah 


where 


(la) 

(lb) 


The  conductivity  and  permittivity  arc  given  by  O  and  e  respectively.  Thu 
permeability  of  all  media  is  taken  to  be  no,  the  permeability  of  free  space. 

The  electromagnetic  noise  from  natural  sources  over  land  is  assumed  to 
produce  a  nearly  vertically  polarized  plane  wave.  A  solution  to  Maxwell's 
equations  for  a  wave  with  this  polarization  propagating  over  a  conducting 
surface  where  the  dirt.tion  of  propagation  is  parallel  to  the  surface  is 
veil  known  and  is  given  by'  J 
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2y  X 

h  -  e 


y  «  0  (2b) 


where 


For  the  frequcnclesi  permittivities,  end  conductivities  considered  in 
this  t.eatraent,  the  quantity  jk/c* «  1.  Under  this  assumption.  Equations 
(3a),  (3b),  and  (3c)  become 


3, 
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Since  k  ia  independent  of  the  conductivity  of  the  guiding  surface,  the 
o 

incident  wagnetic  field  becones 
1 V  4  lkox 

htnc  *  *  y  k  0  (5a) 

"ik2yy  *  ikox 
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vh«r«,  from  Equation*  (4b)  and  (4c) 
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and  1®  the  conductivity  of  the  earth. 

2.3  Contribution  Through  the  Air-Sea  Boundary 

la  this  section  an  approximate  expression  for  the  magnetic  field  in 
the  ocean  due  only  to  the  energy  entering  through  the  air-sea  interface 
v;ill  be  determined. 

Settir  .,  y  equal  to  zero  in  Equations  (5a)  ard  (5b),  it  is  noted  that 

the  incident  magnetic  field  on  the  surface  of  the  e  *th  i*  Independent  of 

the  conductivity  of  the  ear  .  This  is  due  to  the  fact  that  k  is  inde- 

o 

pendent  of  the  conductivity  of  tV  guiding  plane.  Because  of  this,  propa¬ 
gation  in  the  magnetic  fie!  :  can  be  considered  to  be  unperturbed  at  the 
shoreline.  The  magnetic  field  at  the  dp  surface  of  the  ocean  is  then 
given  by 


ikx 

o 


(7) 


To  find  the  field  «’  some  point  in  the  occ  .n,  waves  will  be  considered 

to  ent«r  the  water  normal  to  :he  surface.  This  assumption  Is  based  on 

# 

the  fact  that  propagation  from  a  poor  conductor  to  a  good  conductor  re* 
suits  in  an  angle  of  refraction  which  is  very  nearly  *tro^^«  The  fierd 
it  some  point  In  the  oceen  will  thus  be  of  the  form 


-ikj (y*d) 
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(8) 


where 
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and  is  the  conductivity  of  the  sen  water.  It  is  to  be  noted  that  d, 
the  depth  of  the  ocean,  has  been  introduced  into  the  exponents  in  Equa¬ 
tion  (8).  This  is  done  in  the  interest  of  simplifying  the  final  expres¬ 
sion.  The  subscript  t.  indicates  that  this  is  the  magnetic  field  in  the 
ocean  due  to  the  energy  entering  through  the  Lop  surface  only.  The  quan- 

/g\ 

tity  p,  is  given  by 


-3 

The  conductivities  of  the  earth  and  sea  water  are  given  to  be  10 

(12) 

mhos/meter  and  either  3  or  4  mhos/meter  respectively  •  With  these 
values  is  very  nearly  equal  to  -1.  Continuity  of  h  across  the  bound¬ 
ary  at  y  *  0  requires  that 
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The  final  expression  for  .h  then  becomes 
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2.4  Contribut ion  Th r cvgh  the  Land-Sea  Boundary 

In  this  section  the  rugrctu-  field  in  the  sea  due  only  to  the  energy 
entering  thiough  the  :  I  iv  ocean  boundary  is  determined. 

Consider  an  arbitrary  closed  surface  s'  bounding  a  source  free  region. 
It  can  be  shown  that  a  field  inside  this  region  at  some  point  r  is  given 


*(r)  -  -  T5  f  [c<W--'>!;,  ■<>")  -  Hr’)  |;;,G<r|r*j]  ds*  (10) 

s' 

where  G(r|rf)  is  known  as  the  Green’s  function  which  is  the  solution  to 
the  inhomogeneous  Hemho'tz  equation 

(V2-*-  k2)  G(rjr’)  E  i(te  6  (j  r  -  r'  J);  (11) 

6<i  r  -  r'J )  is  the  Dirac  delta  function,  r*  is  the  coordinate  of  a  point 
source  on  s’,  r  is  the  field  coordinate,  and  the  operator  3/Dn'  denotes 
the  derivative  normal  to  the  surface  s’* 

I  ' 

Referring  to  Figure  2,  it  can  be  seen  that  a  closed  surface  s'  can 
be  selected  such  that  some  of  it's  elements  are  parallel  to  the  z  axis  and 
the  remaining  surface  elements  are  parallel  to  the  x,  y  plane.  The  sur¬ 
face  areas  made  up  of  elements  parallel  to  the  x,y  plane  are  outlined  by 
the  contour  i’.  Setting  'i  ■  h  it  is  noted  that  since  3/dz  ■  0,  the  normal 
derivatives  of  f(r')  and  G(r|r’)  with  respect  to  r'  vanish  on  all  elements 
of  s'  parallel  to  the  x,y  plane.  As  a  result  of  this,  Equation  (10)  re¬ 
duces  to  a  line  integral  along  i',  i.e., 


9. 


h(r>  -  -  |  [g (z  | t  * )  f,,  h(r')  ••  h(r»)  f-~,  G(r)r')  ]  dJt*  (12) 

2  A' 

wheie 


and  x*  and  y'  are  simply  the  rectangular  coordinates  of  the  point  source 
on  the  surface,  s'. 

The  Green's  function  will  be  chosen  such  that  it's  normal  derivative 
will  vanish  along  y  *  -d  for  all  x.  In  addition,  it  is  again  temporarily 
assumed  that  the  ocean  is  a  perfect  conductor  such  that  dh/dn'  vanishes 
for  y  "  -d,  x  *  0.  Under  these  restrictions  the  integral  along  y  *  -d, 

0  *  x  s  +  «  is  zero.  Also,  by  Equations  (3a)  and  (3c),  it  can  be  seen 
that  the  integral  along  the  infinite  arc  is  also  zero  due  to  finite  los- 

i  * 

ses.  Thus  Equation  (10)  reduces  to  an  integral  along  the  y  axi^  from 
y  ■  -d  to  y  -■**.  , 

The  solution  to  Equation  (10)  subject  to  the  restriction  that 
d/dy’  G(rjr')  '  0  at  y  ■  -d  is  well  known  and  is  given  by^^ 


where  H  ^ 
o 

of  the  first 


kind  and 


is  the  zeroth  order  Hankel  function 
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This  Green's  function  represents  the  magnetic  field  at  r  due  to  a 
line  source  at  r'  where  the  line  source  is  described  by  the  Dirac  delta 
function  6(y-y')  times  iu^*  lb's  physical  counterpart  would  be  a  mag¬ 
netic  line  current  source. 

In  order  to  evaluate  Equation  (12)  it  is  necessary  to  know  the  mag¬ 
netic  field  and  its  derivative  with  respect  to  x'  along  y',  the  path  of 
the  integration.  This,  of  course,  is  not  known  exactly.  As  is  often 
the  case  in  problems  of  this  type,  we  shall  approximate  h  over  this  path 
by  assuming  that  it  is  equal  to  the  incident  magnetic  field  which  is  given 
by  Equation  (5b)v  .  Substitution  of  Equation  (5b)  and  Equation  (13)  in¬ 

to  Equation  (11)  yields 


This  integral  cannot  be  readily  evaluated  for  all  values  of  x.  However,  for 

2 

sufficiently  Large  X,  say  x  >  5  6.,  where  6-  *  ’ •y  is  the  skin  depth  in 

the  earth,  it  is  shown  in. the  appendix  that  Equation  (14)  can  be  approxi¬ 
mated  by 


ik^(xtd) 


(15) 


The  assumption  is  now  made  that  waves  enter  the  ocean  in  a  direction  nor¬ 
mal  to  the  bottom  surface.  This,  as  before,  is  justified  by  the  fact  that 


»  O^.  The  field  at  scire  point  in  the  ocean  then  becomes 


The  subscript  b  indicates'  'hat  ‘his  is  the  magnetic  field  in  the  ocean 
due  to  the  energy  entering  through  the  bottom  of  the  ocean  only.  Pfc  is 
the  reflection  coefficient  at  the  top  of  the  ocean  and  is  given  by 


This  reflection  coefficient  is  also  very  nearly  equal  to  -1.  Requiring  h 
to  be  continuous  across  the  boundary  yields 


s. 


a  +  Rta 


The  final  expression  for  now  become * 
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N. 


no 


This  derivation  has  neglected  the  vertical  portion  of  the  earth-ocean 


interface.  This  la  justified  by  the  fact  that  the  skin  depth  in  the  ocean 
turns  out  to  be  always  considerably  less  than  the  shortest  distance  from 


12. 


shore  for  t uioh  the  final  expression  is  vjlld. 

2.5  Total  Magnetic  Field  in  the  Oce.in 

The  total  magnetic  field  In  the  ocean  is  givei  by  fhe  sum  of  Equa- 


tion  (9)  and  Equation  (16)  whh.h  gi'es 


total 
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■t-  h.  (17) 
d  bs 


Setting  pfc  and  equal  to  -1  allows  Equation  (17)  to  be  reduced  to 


Sinkj  (y+d)  Sinkjy 

htotal  *  Si^a  hts  "  sln^jd  hbs 


where  h  and  h.  are  given  by  Equation  (?)  and  Equat  ion  0-*)  respectively. 

ts  b* 
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III.  THE  ELECTRIC  FIELD 


The  vertical  and  horizontal  components  of  the  electric  field  are 
found  in  the  first  part.  of  this  chapter.  The  two  field  components  are 
then  compared  with  respect  to  their  relative  magnitudes. 

3.1  The  Electric  Field  Components 

Now  that  an  expression  for  h  hot*  been  obtained)  the  electric  field 
is  easily  obtained  by  substituting  Equation  (13)  into  Equations  (la)  and 
(lb).  The  horizontal  component  of  the  electric  field  thus  becomes 


k.  p 
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COskjCy+d) 


Sinkjd 


coakjy 

Sinkjd 


ts  Sink 
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and  the  vertical  component  becomes 


(19) 


Sinkjd  I  »x  (JO) 

where,  again,  usa  has  been  made  of  the  feet  that  mt./*j  ^  1* 

figaMXittS 

Although  tho  vortical  and  horizontal  components  of  tho  electric  field 
have  been  determined,  the  dlreotlon  of  tho  not  eleetrie  vector  la  not  im- 

mediately  apparent.  In  this  section  the  magnitude*  of  the  two  components 

\ 

of  tha  electric  vector  will  be  compared. 

Consider  first  the  problem  of  the  infinitely  deep  ocean.  Due  to  the 
absence  of  the  sarth-ocaan  interface,  there  would  be  no  energy  radiated 
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in  the  +y  direction, 
with  k^y  replaced  by 


in.*  solution  would  then  reduce  to  Equation  (5b) 
k,.  I  he  magnetic  field  would  then  become 
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where,  as  before,  th*  ;  •  v;  ibution  through  the  vertical  earth-ocean 
interface  has  been  ig-.;rei.  From  Equa-  ions  (la)  and  (lb),  the  electric 
field  components  now  become 
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(22b) 


Dividing  Equation  (2/.b)  by  (22a)  and  taking  the  absolute  magnitude  yields 
E  k  ft wT 

1  ^  i  ■  i  *7 1  -  ysf 

-12 

Setting  «Q  *  8.85  x  10  farada/metcr  and  ffj  •  4  mhos/meter,  Equation 
(23)  bacomea 

|  |  »  l.S  fm  x  lo”® 

V 

It  la  apparent,  therefore,  that  for  frequencies  of  10*°  c/s  and,. less, 
the  vertical  component  of  the  electric  field  may  be  neglected  in  compart* 
•on  to  its  horltontal 

This  roeult  is  no  surprise  since  it  already  has  been  pointed  out 


that  tha  waves  are  entering  the  oceen  with  a  direction  of  propagation 
normal  to  the  surface.  1*.  is  not  so  obvious,  however,  that  this  result 


o. 


should  hold  for  an  c.e-in  with  a  finite  depth  for  it  might  be  argued  that, 
at  some  point,  the  upward  traveling  waves  and  the  downward  traveling 
waves  might  have  x  components  of  their  respective  electric  fields  which 
combine  destructively  wh-reas  their  y  components  combine  constructively. 

It  shall  nov;  be  ••■how-.  thj‘-  the  vertical  component  can  still  be  ne¬ 
glected  under  the  re -r rlctior-s  of  the  model  used  in  this  study. 

Since  no  values  cf  d  less  than  5  ocean  skin  depths  are  used,  the 
approximation  is  made  that  reflect  ions  at  the  earth-ocean  and  air-ocean 
interfaces  may  be  neglected.  The  magnetic  field  may  then  be  written  as 


-ik.y  ik  (ytd) 
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(24) 


Substitution  of  Equations  (7)  and  (15)  in-o  Equation  (24)  yields 


h  mt  e 
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Use  of  Equations  (la)  and  (lb)  then  give  ttus  electric  field  as 
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where,  aino  k  «  k„,  k,  has  been  neglected.  Taking  the  absolute 
0  z  < 


I 


i 


magnitude  of  the  ratic  cf  E  to  E  results  in 

y  * 


i  ♦ 


ik^(xtd)  *  ikj(2y+d)~  ikQX 

J  k  n(  xtd)TIO  2y+d  )  -”HTx 
2  l  o 


(26) 


where  use  has  been  madr  of  the  fact  that 


1  1 
Defining 


Equation  (26)  is  rewritten  as 


S  -Ae 


ie 


1*  Ae 


i© 


where 


A  - 


|2"k?x 


.  •  -  J  ♦  <1M>  |  J*  |  (*»«•«♦»>  |^i 


k 


(2>>d) 


♦  1 


(fcl 


hr? 

h  I 

V’! 


'> 


<2y*d)  1*  «♦»» 


(27) 


(2S«) 


<28b) 
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S  *  k0/k2 


(28c) 


Equation  (27)  now  becot.u  a 


S  -  A  f  Cos  (a> lg)  +  lStt-^lf1)} 
1  +  A  {  Co?.  (O'i?)  -r  iSin(a+tP)J 


(29) 


Using  the  identities 

Cos(fffi 0)  ■  CosaCoshP  -  i5inSSinhp 
Sin(OreiP)  *  SinOCoshP  ♦  iCosOSinhP 

Equation  (28)  can  be  rewritten  as 


M 


{SR  -  ACosflr  [CoahB-SinhS] }  ♦!  {Sj+ASinofSinhP-CoshB] J 
(1  "<♦  ACostt  [CoahS-SinhPj}  -i  t  ASinfl^SinhP-CoshPj  J 


where  Sg  ♦  iSj.  •  S 

Equation  (30)  can  be  further  simplified  by  making  the  substitution  . 

♦  [  CoshP  -  SinhP  ]  •  ♦  e’^ 

Then  Equation  (30)  finally  becomes 


fSR  -  ACosbre*®}  ♦!  (Sj  -AS»n«c'*} 

— 1  *  ■  1 —■■■■<■ . . 

{la  ACos«e  }  ♦!  (ASlne”®} 


(31) 


In  order  to  reduce  thv  corplcxlty  of  Equation  (31)  the  restriction 
that  I  *  Q  will  be  Imposed  a-d  the  consequence*  of  this  restriction 
will  later  be  examined.  8  Is  defined  by  Equation  (28b). 


Setting  P  -  0  and  '.o'.r.g  that  lor  the  parameters  used  in  this  model 
A  s  0.15,  the  denominator  of  Eqearion  (3i)  can  now  be  represented  in  the 
Inequality 


{  1  +  ACc-,Ci  }  M  (  Asi: a  }  >  (1-0.15)  -  0»85 


Under  the  same  rt-stri.  *  *.h..  numerator  has  an  upper  bound  given  by 


t  Sjj-ACosa}  +i  l  S,- AS>nCr}j  *  A(Cosor  +i  Sindr)  S  JT  (0.15)  -  0.212 


where,  since  S  “  SR  *  5.2  x  10  ,  they  have  been  neglected.  Under  the 

restriction  then  tha*  I  ■  0.  ’he  maximum  value  that  |B  /B  I  can  take  on 

1  y  x1 

is  given  by 


^l*MUI‘MKIK  W 


Now  since  A  *  0.15,  i*  is  ob.irus  from  Equation  (31)  that  Equation  (32) 
holds  also  for  0  >  0.  ^  »  4  mhos/mcter  and  •  10*^  mhos/meter 

ths  inequality  (32)  bcr.omes 


.  j  ^  |  *  4  x  10~3 


In  order  to  gcairc*.?  St  si  ly  dvfir.c  the  ugion  in  which  this  inequality 
holds.  Equation  (23b)  is  *M\\d  f-r  y  with  P  •  0.  This  results  in 


y  *  *  i  -  i  |  |  *  *  2  *  1  <**d> 


19. 


yields  a  minimum  of  -£ .42  fer  B«  He-  it  Is  now  only  necessary  to  con 
aider  8  in  the  range  0  *  8  >  -.V.'»2, 

ft 

As  B  ranges  from  0  to  e  ranges  frero  1  to  11.3  and  cos« 

ranges  fiotr.  0,95  *.o  -  0-i  ij  t" d  passes  through  zero  only  once.  The 
lower  bound  of  the  de'emr.  a*.,r  of  Ecrvir-n  (31)  is  then  giver,  by 

j  {  1  +  ACcsOfc*'®  }  j  *  1  +  A  <U,3)  (-0.276) 
but  A  £  0.15,  then 


1  +  ACosore”^  J 


Si-  0.466  *  0.534 


The  numerator  of  Equation  (31)  has  an  upper  bound  given  by 


j  ACosOfe"8-iASir0ft'®  |  s  (0.15)  (11.3) 

where,  as  before,  SR  and  Sj  have  bc-en  neglected.  It  is  concluded  then, 
that  when  0  >  M  -2.42, 


The  restrictions  or  these  results,  that  is  that  outside  the  cross 
hatched  region  |Ey/Ex (  s  4  x  10  3  and  inside  the  cross  hatched  region 
lEj/Exl  S  7*2  x  10  ,  are  first,  that  the  ocean  depth  be  at  least  56, 
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ii 

.1 


so  as  to  validate  the  assumption  that;  reflections  at  the  earth-ocean  and 
air-ocean  interfaces  ar\  negligible,  and  secondly  that 

0.097  <  A  s  0.1: 

This  is  equivalent  tc  'vqutr:ln§  that 

562  <  x  «  1262 

For  x  less  than  bfi^s  the  far  field  approximation  given  by  Equation  (15) 
no  longer  holds  and  for  x  greater  than  1262»  the  contribution  through 
the  earth-ocean  interface  is  eliminated  to  such  an  extent  that  the 
profiles  approach  a  simple  exponential  form. 

The  variation  of  A  with  a is  extremely  slow,  such  that  a  change  in 
0£  of  plus  or  minus  an  order  of  magnitude  has  a  negligible  effect  on 
these  results. 

Comparing  the  upper  brunds  on  js  /E  |  it  is  noted  that  inside  the 
cross  hatched  region  in  Figure  3,  the  upper  bound  is  a  factor  of  18  larger 

than  the  upper  bound  outside  of  the  cross  hatched  region.  This  points  out 

» 

the  possibility  that  very  near  the  shore,  in  the  region  where  this  model 
is  no  longer  valid,  |Ey|  may  be  of  an  appreciable  size  with  respect  to 

|sx|. 


IV.  NUMERICAL  RESULTS 


Numerical  values  for  the  magnitude  of  the  x  and  y  components  of  the 
electric  field  were  determined  by  digital  evaluation  of  Equation  (19) . 

The  field  is  normalized  by  multiplying  by  The  results  are  plotted 

in  Figures  4  to  16  as  a  function  of  depth  and  distance  from  shore  for  fre¬ 
quencies  of  10v  60,  100 1  and  1000  oycles/sec.  Depths  range  from  50  meters 
to  800  meters  and  distances  from  shore  range  from  3  km  to  50  km.  The  depth 
into  the  ocean  is  plotted  as  the  absolssa  and  the  various  values  of  x  (dis¬ 
tance  from  shore)  determine  a  family  of  ourves. 

In  Figure  17  a  comparison  is  made  with  some  experimental  points  de¬ 
termined  at  a  frequency  of  60  cycles  /seo.  and  a  depth  of  J00  meters 

The  value  of  x  has  a  lower  limit  which  is  determined  by  the  approxi¬ 
mation  used  to  evaluate  the  integral  in  Equation  (l4).  This  is  essen¬ 
tially  a  far-field  approximation  which  requires  x  to  be  no  smaller  than 

about  56_  where  6_  -  \l  — - —  is  the  skin  depth  of  the  earth.  6, .  the  skin 

2  2  V  MUgCg  t  1 

depth  in  the  ooean,  and  6^  are  plotted  for  convenience  as  a  funotlon  of 

♦ 

frequency  in  Figures  19  and  18  respectively. 

The  contribution  to  the  field  due  tc  the  energy  passing  through  the 

air-ocean  boundary  only  is  shown  as  a  dotted  line.  This  enables  one  to 

■ 

see  store  explicitly  the  effeot  that  the  energy  passing  through  the  horl- 
sontal  land-ooean  interface  has  on  the  total  field. 


V.  SUMMARY  AND  CONCLUSIONS 


5 


The  electromagnetic  field  in  the  ocean  resulting  from  excitation  by 
natural  sources  over  land  has  been  investigated  with  the  intent  of  gaining 
some  insight  as  to  the  mechanism  by  which  the  energy  actually  enters  the 
water.  The  horizontal  and  vertical  components  of  the  electric  field  are 
plotted  for  several  values  of  the  parameters  involved  and,  in  one  case, 
a  comparison  is  made  with  some  experimental  measurements. 

Before  drawing  any  conclusions  it  should  be  stated  that  both,  the 
land-sea  model  and  the  configuration  of  the  incident  field,  which  have 
been  chosen  are  only  a  very  rough  approximation  of  the  real  physical 
situation.  Nevertheless,  the  approximations  used  are  felt  to  be  justi¬ 
fied  by  both,  the  complex  nature  of  the  problem  and  the  reasonable  agree¬ 
ment  with  some  experimental  results. 

Upon  examination  of  the  curves,  it  is  seen  that  the  contribution  to 

the  net  field  which'  enters  through  the  land-ocean  interface  is«quite 

prominent  in  the  region  near  the  shoreline.  As  the  distance  from  shore 

increases,  this  contribution  becomes  less  and  lass  noticeable.  This  is 

due  to  the  fact  that  kx,  the  wave  number  which  governs  the  propagation 

in  the  x  direction  for  the  contribution  through  an  air-ocean  interface,  ' 

has  a  small  Imaginary  part,  whereas  the  attenuation  of  the  contribution 

0  ‘ 

through  the  land-ocean  interface  is  quite  high  due  to  the  conductivity 
of  the  earth. 

* 

The  affect  of  the  depth  of  the  ocean  on  the  electric  field  le  also 
q(uita  noticeable.  For  the  region  near  the  shoreline  the  curves  show  the 


bottom  contribution  more  distinctly  as  the  depth  of  the  ocean  increases. 
This  is  understandable  since  for  a  deep  ocean  the  energy  entering  the 
top  surface  will  be  greatly  attenuated  before  reaching  the  bottom.  The 
contribution  through  the  b:ttcm  will  therefore  be  more  in  evidence* 

The  comparison  with  experiment  made  in  Figure  17  is  rather  hampered 
by  a  lack  of  data  points.  It  is  to  be  noted,  however,  that  an  Increase 
in  the  field  intensity  13  clearly  evident  in  the  experimental  results. 

The  discrepancy  between  the  experimental  points  and  the  theoretical 
curves  could  be  due  to  any  or  all  of  the  following:  The  crudeness  of 
the  mathematical  model,  differences  between  actual  and  used  valuta  for 
conductivity,  conductivity  gradients  in  the  earth,  and  reflections  from 
boundaries  in  the  earth  below  the  ocean  bottom. 

In  future  work  a  more  complicated  model  is  to  be  considered.  Mora 
specifically,  the  earth  beneath  the  ocean  will  ba  assumed  to  be  composed 
of  two  layers.  It  la  hoped  that  the  model  will  give  even  better  agree¬ 
ment  with  experimental  results.  Considering  the  problem  in  reverse,  this 
double-layer  model  could  produce  a  method  for  determining  not  only  the 
conductivity  of  the  earth  under  the  ocean,  but  also  the  thickness  of  the 
top  layor. 
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n - 2  X  3°2 

J x*(y,+d)*  *  x 

.  r  4  2‘2 

Since  the  va !■.:;?  cf  d  considered  are  much  less  than  its  presence 
in  the  argument,  ccei  net  afftnr  :l>i  approximation  to  any  great  extent. 

With  these  approximations,  the  Hankel  functions  are  no  longer  func¬ 
tions  of  y*  and  therefore  may  be  mo*ed  outside  of  the  integral  sign 
leaving  a  simple  exponential  to  integrate.  Upon  integrating,  Equation 
(35)  becomes 


lk2d  r  1 

-^7  [kA(1)  <V> 4  lk2  Bia)  V’J 


(36) 


(I) 


Factoring  ikjHj  * (k^x)  outaide  the  brackets  reaulta  in 


h  1  ‘V  .I  tf),k  . 

''b*  *  -  2  •  h  '“g*' 


k 

M  t 


H^(1)(k,x) 


x  o  X- 
2  H, '"'(k,x) 


■or 


(37) 


given  at 


The  large  argument  approximations  for  the  Hsnkel  functions  are 

(4) 


I 


Dividing  Equation  (33a)  by  Equation  (3?b)  yields 


H  (1)(K,x)  i  | 
o  2  ...  2  _ 


— 7H — - —  &  e  *  1  for  large  k_x 

H1(1,<V)  2 


This  considerably  simplifies  Equation  (32)  i.e.. 


Noting  that  «  k2>  Equation  (39)  finally  become 
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